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ABSTRACT 

We use the Kossakowski-Lindblad-Davies formalism to consider an open system de¬ 
fined as the Markovian extension of one-mode quantum oscillator S , which is perturbed 
by a piecewise stationary harmonic interaction with a chain of oscillators C. The long¬ 
time asymptotic behaviour of various subsystems of S + C are obtained in the framework 
of the dual IU*-dynamical system approach. 
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1 Introduction 

A quantum Hamiltonian system with time-dependent repeated harmonic interaction was 
proposed and investigated in TZl . The corresponding open system can be defined through 
the Kossakowski-Lindblad-Davies dissipative extension of the Hamiltonian dynamics. In 
our previous paper jTZlj the existence and uniqueness of the evolution map for density 
matrices of the open system are established and its dual kP*-dynamics on the CCR C*- 
algebra was described explicitly. 

The aim of this paper is to apply the formalism developed in jTZlj to analysis of 
dynamics of subsystems, including their long-time asymptotic behaviour and correlations. 

Let a and a* be the annihilation and the creation operators defined in the Fock space 
& generated by a cyclic vector Q (vacuum). That is, the Hilbert space & is the com¬ 
pletion of the algebraic span J^ fin of vectors {(a*) m Q} m ^ 0 and a, a* satisfy the Canonical 
Commutation Relations (CCR) 

[a, a*] = 1, [a, a] = 0, [a*, a*} = 0 on ^ fin . (1.1) 

We denote by {J^k} k =o the copies of & for an arbitrary but finite JV6N and by 
the Hilbert space tensor product of these copies: 

N 

jfW ■=(g)j?’ k = N +d . (1.2) 

k =o 

In this space we define for k — 0,1, 2,..., N the operators 

fe fc :=l®...<g)l<g)a<g)l<g)...<g)l, b* k := 1(8). ..<8)1(8) a* (8)1(8). ..<8)1, (1.3) 

where operator a (respectively a*) is the (k + l)th factor in (II.3|h They satisfy the CCR: 

[h,b* k ,] = 6 ktk ' 1, [b k ,b k >} = [b* k ,b* k ,} = 0 (k, k' = 0,1,2,..., N) (1.4) 


2 




















on the algebraic tensor product (^ r &n }^ N+1 ' > ■ 

Recall that non-autonomous system with Hamiltonian for time-dependent repeated 
harmonic perturbation proposed in [ TZ] has the form 

N N 

H N (t ) := Eb* 0 b 0 + e ^ b* k b k + 77 ^ X[(fc-i)r,fcT)(<) (&S&* + &fc&o) • (1-5) 

k =1 /c=l 

Here t G [0, iVr), the parameters: r,E,e,r) are positive , and X[x,</)(•) is the characteristic 
function of the semi-open interval [x,y) C R. It is obvious that H^(t) is a self-adjoint 
operator with time-independent domain 

N 

Vq = P| dom ( b* k b k ) C . (1.6) 

k =o 

The model (11.51) presents the system S + Cn, where S is the quantum one-mode cavity , 
which is repeatedly perturbed by a time-equidistant chain of subsystem: Cn = Si + S 2 + 
... + Sn- Here {<Sfc}fc>i can be considered as “atoms” with harmonic internal degrees 
of freedom. This interpretation is motivated by certain physical models known as the 
“one-atom maser” [ BJM j. [NVZj . The Hilbert space J^s •%> corresponds to subsystem 
S and the Hilbert space to subsystems S k (k — 1,..., N), respectively. Then (II .2[) is 

N 

Jd’W=je’s®'%c N , Mc n -.= { (1.7) 

k= 1 

By (1 1.5 p only one subsystem S n interacts with S for t G [{n — l)r, nr). In this sense, the 
interaction is tuned m The system 5 +Cat is autonomous on each interval \{n — l)r, nr) 
governed by the self-adjoint Hamiltonian 

N 

H n \= E b* 0 b 0 + Kh + V ( K b n + K b o) , n = 1, 2,..., N , (1.8) 

k =1 

on domain V 0 . Note that if 

rj 2 ^ E e , (1.9) 

Hamiltonians (11.51) and (II. 8ft are semi-bounded from below. 

We denote by (J^^) the Banach space of the trace-class operators on . Its 

dual space is isometrically isomorph to the Banach space of bounded operators on : 

~ The corresponding dual pair is dehned by the bilinear functional 

(0| A)jg, w = Tr jp(N) (0 A) for (0, A) G x C(^ N) ) . (1.10) 

The positive operators p G € 1 (Jf p '- N l) with unit trace is the set of density matrices. 
Recall that the state u> p over is normal if there is a density matrix p such that 

^p(-) = (p\ • (1-11) 
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1.1 Master equation 

To make the system S + open, we couple it to the boson reservoir 1Z, [AJP3]. More 
precisely, we follow the scheme (S + 1Z) + Cn, he. we study repeated perturbation of the 
open system S + TZ [ NVZj . 

Evolution of normal states of the open system (S + TZ) + Cn can be described by the 
Kossakowski-Lindblad-Davies dissipative extension of the Hamiltonian dynamics to the 
Markovian dynamics with the time-dependent generator [Ah] , [AJP2| 

La(t)(p) := -i[H N (t),p] + (1.12) 

+c(p)-1(s*(i)p+pc*(i)), 

for t ^ 0 and p <G dom L a {t) C <£i(Ji Here the Erst operator Q : p (->• Q(p) G 
in the dissipative part of (11.12)1 has the form: 

Q(-) = a- b 0 (•) b* 0 + a+ b* 0 (•) b 0 , a T ^ 0 , (1.13) 

and the operator Q* is its dual via relation (Q(p) \ A)j^( N ) = (p 

Q*(-) = a-b* 0 (-)bo + <J + b 0 (-)b* 0 . (1.14) 

By virtue of (jl .5)1 . for t G [(n — l)r, nr), the generator (II, 12)1 takes the form 

L„(p) :=-i{H„p}+ a(p)-i(a*(l)p + pC*(l)). (1.15) 

The mathematical problem concerning the open quantum system is to solve the Cauchy 
problem for the non-autonomous quantum Master Equation |lA.IP2) 

dt pif) = L a (t)(p{t )) , p(0) = p. (1.16) 

For the tuned repeated perturbation, this solution is a strongly continuous family {T ° 0 } t > 0 , 
which is defined by composition of the one-step evolution semigroups: 

rn(T _ nr\(j rri(j rji(j rri(j 

1 t,0 — 1 t,(n-l)r 1 n-l * * * 1 2 1 1 ? 

where t = (n — l)r + v(t), n ^ N, v(t) < r. Here we put 

7T:=r?M, TZ(s):=e- L -* (O 0), (1.17) 

and then = T°(y(t)) holds. The evolution map is connected to solution of the 

Cauchy problem (11.161) by 

T t,o : P ^ P(t) = T t,o(p)- (1-18) 

The construction of unique positivity- and trace-preserving dynamical semigroup on 
for unbounded generator (jl , 15)1 is a nontrivial problem. It is done in jTZlj 
under the conditions (11.91) and 

0 ^ cr + < cr_ . (1.19) 

for the coefficients in (11.131 II.14ft . Then, {T% (s)} s ^o for each k (11.171) is the Markov 
dynamical semigroup, and (11.181) is automorphism on the set of density matrices. 
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1.2 Evolution in the dual space 

In order to control the evolution of normal states, it is usual to consider the IW-dynamical 
system {Tf 0 *} t ^ 0 ), where are weak*-continuous evolution maps on 

the von Neumann algebra ~ |AJP1] . They are dual to the evolution 

(ll.ISp on Ciby the relation (ll.lOlh 

(TlM\A) jem = {p\Tl‘(A)) jem for (p,A) e x C{.X’ {N) ) , (1.20) 

which uniquely defines the map A t—>■ T t a 0 *(A) for A e The corresponding dual 

time-dependent generator is formally given by 

L:(t)(-)=t[H N (t),-}+ (1.21) 

+ a*(')-l(a*(i)(-) + (-)c*(i)) for t>o. 

When t G [(k — l)r, hr), the above generator has the form 

i;y.) = i[H„, •] + e*(-)-l(e*(i)(o + (oe*(i)). (1.22) 

We adopt the notations 

TV = TW . r,y„_ 1)r = r„>(f))*. and irw* ~(■>«), (1.23) 

dual to (11.17|) for t = (n — 1 )r + z/(f), n ^ N, v(t) < r. Then, we obtain 

r,7(d) = Tf *t 2 ct *... r;*i); ( ; 1)T (A) for a e c {.^ (JV) ). (1.24) 


Let (or CCR(C) ) denote the Weyl CCR-algebra on &. This unital C'*-algebra 

is generated as operator-norm completion of the linear span si w of the set of Weyl oper¬ 
ators 

w(a) = e i ^ a) (ae C), (1.25) 

where <E>(a) = (aa + aa*)/V% is the self-adjoint Segal operator in &. [The closure of the 
sum is understood.] Then CCR (II. lj) take the Weyl form 

w(ai)w(a2) = e ~ lIm ( aia2 )/ 2 uj{ai + a 2 ) for a 1 ,a 2 G C . (1.26) 


We note that is contained in the C*-algebra £(lF) of all bounded operators on &. 

Similarly we define the Weyl CCR-algebra si C over This 

algebra is generated by operators 


w(C) 



for ( 


(( o\ 

Cl 


e c N+1 . 


\C n) 


(1.27) 
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By (11.31) . the Weyl operators (II ,27j) can be rewritten as 


W(£) = exp[i((C,6) + (bX))/V2\ , 

where the sesquilinear form notations 

N N 

<C, &> := £ CA, (6. C):=£CA* 


(1.28) 


(1.29) 


J=0 


J=0 


are used. Let us recall that is weakly dense in [AJPlj . 

Explicit formulae for evolution operators (II.23jl acting on the Weyl operators has been 
established in |TZlj . For n — 1,2,... .N, let J n and X n be ( N + 1) x (N + 1) Hermitian 
matrices: 


( Jn)jk ~ 


1 (j — k — 0 or j — k — n) 
0 otherwise 

/ 


(1.30) 


\E-e )/2 (j, k) = (0, 0) 

— (E — e)/2 (j, k) — (n, n) 

(■ x n)j k = V (j, k) = (0, n) ■ 

V Cb k) = (n, 0) 

0 otherwise 


(1.31) 


We dehne the matrices 


Y n :=el + 


E-e 


Jn + X n (n — 1,.. •, N ), 


(1.32) 


where / is the (N + 1) x (N + 1) identity matrix. Then Hamiltonian (II. 81) takes the form 

N 

H n — ( Y n)ik,b*b k . (1.33) 


j,k =0 


We also need the (N + 1) x (N + 1) matrix P 0 defined by (P 0 )jfe = $jo$ko (. j, k = 
0,1,2, ...,1V). Then one obtains the following proposition which is proved in [ TZlj : 


Proposition 1.1 Let n = 1, 2,..., N and £ 6 C^ -1-1 . Then for s ^ 0, the dual Markov 
dynamical semigroup HI. 25 1) on the lTe?/t C*-algebra has the form 


where 


Tr(s)(w(0) = KAOw(K(s)() , 


^L(C) := exp 


1 cr_ + cr + 
4 cr_ — cr_)_ 


«c,0-WWCCCWO) 


(1.34) 

(1.35) 
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and 


K( s ) = ex P 


is[Y n + i 


. a- — a+ 


Pr 


(1.36) 


under the conditions m and HI. 19) . Therefore, the k-step evolution (t = kr , k ^ N in 

«C,C>- {Uf...uzc uz...uiQ) 


( [ 1.24\ ) ) of the Weyl operator is given by 

n T i(W(Q) = exp[ - g - + 

L 4(cr_ — cr_|_) 

x W(U1...VIQ, 

where ‘ = Tf *T, * ... T %' and U’ := U;(r). 


(1.37) 


Remark 1.2 27ie explicit expression of the matrix U%(t) in H1.36\) is given by Uf{t) = 
e lte V°(t), where 


g a (t)z a (t) 8 k o + g rT (t)w (T (t) 8 kn 
( v n(t))jk = { g ff (t)w a (t) 8 k o + g a (t)z a (—t) 8 kn 

8jk 


(J = 0) 

(j = n ) 

( otherwise ) 


(1.38) 


fTere E a := E + i (cr_ — a + )/2 and 
g a (t) := e it{E °~ e)/2 , w a (t ) := 


2^ 


z CT (f) := cos t\/ A ^ + rf + 


V(£ ff -e) 2 + 4 t7 2 
i(E a - e) 


sin t 


sin t 


(E 0 


+ 7/ 2 


- e) 2 , „ 2 


+ 77 - 


(1.39) 


(1.40) 


VIB.-eP + V 

IVoie that the relation z a (t)z a (—t) — w a {t ) 2 = 1 holds for any a± ^ 0, whereas one has 
|( 7 <T (t)| 2 (|£ <T (t)| 2 + \w a (t)\ 2 ) < 1 and z a (—t ) ^ z a (t) for 0 ^ cr+ < cr_. 


Hereafter, together with \1.3 7| ) we also use the following short-hand notations: 

9 a = g a (T), w a :=w a (T), z a = z°(t) and V" := V°(r). (1.41) 


Remark 1.3 Dual dynamical semigroups (1.34 ) an d the evolution operator ( 1.37 ) are 
examples of the quasi-free maps on the Weyl C*-algebra. Using the arguments of IDVV f. 
we have shown in !TZ1\J that they can be extended to the unity-preserving completely 
positive linear maps on under the conditions 111.91) and 1)1.19) . 


The aim of the rest of the paper is to study evolution of the reduced density matrices 
for subsystems of the total system (S + TV) + C]y. 

In Section [2j we consider the subsystem S. This includes analysis of convergence to 
stationary states in the infinite-time limit N —y 00 . We also perform a similar analysis for 
the subsystems S + S m and S m + S n . Section [3] is devoted to a more complicated problem 
of evolution of reduced density matrices for finite subsystems, which include S and a part 
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of Cat. This allows us to detect an asymptotic behaviour of the quantum correlations 
between S and a part of Cat caused by repeated perturbation and dissipation for large N 
in terms of those for small N with the stable initial state. 

For the brevity, we hereafter supress the dependence on N of the Hilbert space M^ 
as well as of the Hamiltonian H N (t) and the subsystem Cm, when it will not cause any 
confusion. 


2 Time Evolution of Subsystems I 


2.1 Subsystem S 

We start by analysis of the simplest subsystem S. Let the initial state of the total system 
S + C be defined by a density matrix p £ £i(Ms <E> Mz). Then for any t ^ 0, the evolved 
state 4(') on the Weyl C'*-algebra stf(Ms) °f subsystem S is given by the partial trace: 

4(4 =Vp(t)(A® l) = Tr J r se) je l;N (T t * 0 (ps®pc)A® 1) for 16^(4), (2-1) 

where p(t) = T£ 0 p and 1 £ srf(Mc). Recall that for a density matrix p £ (Ms ® Mz), 
the partial trace of q with respect to the Hilbert space Mz is a bounded linear map 
Tr^. : Q H > Q 6 £i(Ms) characterised by the identity 

Trjr s ®J%(e( A ® !)) = Tr^ 5 (pH) for A £ C(Ms) ■ (2.2) 

If one puts 

p s (t) :=Tr,^ c (T t a 0 (p)), (2.3) 

then one gets the identity 

(4(H) = Tr je s (p s (t) A ) =: w ps(t) (H) , (2.4) 

by (12.11) . i.e., ps(t) is the density matrix defining the normal state 4- 

In the followings, we mainly consider the initial density matrices of the form: 

N 

p = Ps® Pc for ps = Po , Pc = (§Qpk with Pl = p 2 = ... = p N . (2.5) 

k =i 


Note that the characteristic function E us : C —* C of the state ojg on the algebra (Ms) is 


E u S W = us(w(8)) 


( 2 , 6 ) 


and that (12.6(1 can uniquely determine the state ujg by the Araki-Segal theorem (AJPlj . 
Lemma 2.1 Let A = w(9). Then evolution of H2.1\) on the interval [0,r) yields 

\9\ 2 cr + cr + 


E^ s (0) = exp 


4 <j— — <J-j_ 


(l-Kt^wr-i^Kwr) 


> * e [o, t) . 


(2,7) 










Proof : By (11.271) . we obtain that W(9e ) = w{9) <E> 1 <E> • • • <E> 1 for the vector e = 
*(1, 0,..., 0) G c Ar+1 , where *(...) means the vector-transposition , cf (jl.27[) . Then (12.lft - 
(12.4p yield 

t4(tu(0)) = u p(t)(w(0) ® 1 ® ® 1) = ^p 5 w(^(^)) • (2-8) 

By virtue of duality (11,20j) and (11.371) for k = 1, we obtain 
"P5(*) W)) = ^ps^pc ((T;;w)(ee)) = 

= exp f- l^£z-t£±(l - <t7fWe,t7f(«)e))lu l8 » (W(St^(i)e)) . 

“t (T — (J _j_ J 3 

Taking into account (II,38jl and (12.61) . one obtains for (j2.8[) the expression which coincides 
with assertion (12.71) . □ 

Similarly, for t = mr we obtain the characteristic function 

W 2a - + a+ (l-(U?...U°e,U;.--U°e)) 

(2.9) 




W="fls®«:(%W0e)) = exp 


4 cr_ — cr + 

x (tr(»f/r..vM) = 


= exp 


I /O 12 | 2 V 

4 ~ {u ’■ ■ ■ Ke ’ u °- flv R®w• • • u Wj», 


t=o 


where we have used (ll.27j) and (11.371) . By (jl.38[) we obtain 

{ e lArT£ (/(r)/(r)) m (fc = 0) 

e lNTe g a {j)w a { T )z a { T )) m ~ k (1 ^ k ^ m) (2.10) 

0 (m < k ^ N). 

Then taking into account \g a z a \ < 1 (Remark 11.21) . we find 

(e,e)-{U;...U^e,Ur...U^e) (2.11) 

By setting m = N, (12.01) . (I2.9j) - (12. 1 ip yield the following result. 

Lemma 2.2 The state of the subsystem S after N-step evolution has the characteristic 
function 

E u^(0) = uj ps ( Nt) (w(Q)) ( 2 . 12 ) 


= exp 


|0| 2 cr_ + cr + 
4 cr_ — cr + 


(i-i^n i 


,cr,cr 12 


\g a w 


1 - \g a z a \ 2 


N 

■xu K (w( e iN '’(g’) N (z°) N e)) J] u n (vKe iN ™(g’’) N - t+1 (z’’) N - k w’0)) ■ 


k =1 
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To consider the asymptotic behaviour of the state ujff T for large N, we assume that 
the state co Pk on g/(L?) is gauge-invariant , i.e., 

e ~i4> a * a p k e i*a*a = pk ^ e ( 2 . 13 ) 

for each component of the initial density matrix pc (12,5ji . 

Theorem 2.3 Let u> Pk be gauge-invariant for k — 1, 2,..., N and suppose that the product 

OO 

D(e):=Y[uj n (w((g’zye)), ( 2 . 14 ) 

s=0 

converges for any 9 e C and let the map R 3 r H > D(r 9) e C be continuous. Then for 
any initial normal state cAj(-) = co Po (-) of the subsystem S, the following properties hold, 

(a) The pointwise limit of the characteristic functions 112.12\) exists 


E*(9) = Jim uj ps ( Nt )(w(9)) , 9 e C . 


(2.15) 


(b) There exists a unique density matrix pf such that the limit \2.15\ ) is a characteristic 
function of the gauge-invariant normal state: E*(9) = uj p s(w(9)). 

(c) The states {adT r } rn ^i converge to uj p s for m —> oo in the weak*-topology. 

Proof: fal Bv (|1.25p and bv the gauge-invariance (|2.13p. one gets u> n (w(e l< ^9)) = u Pk (w(9)) 
for every (f G R. Hence, for 1 ^ k ^ N the characteristic functions E Upk ( 9 ) depend only 
on |0|, and we can skip the factor e lNre in the arguments of the factors in the right-hand 
side of (12.121) . Note that for IV —> oo the factor u j po converges to one, since the normal 
states are regular and \g (T z a \ < 1 (see Remark 11.21) . Hence, the pointwise limit (12.1511 
follows from (j2.12jl and the hypothesis (I2.14p . It does not depend on the initial state u po 
of the subsystem S and the explicit expression of (12.1511 is given by 


E*{9) = exp 


\9\-o + a + 

b 

9 

b 

2 


4 a- — a + 

1 - \g a z a 

V J 


D(g ff w ff 9). 


(2.16) 


(b) The limit (12.16(1 inherits the properties of characteristic functions E u rnr (A)=wrW)): 
(z) normalisation: i?*(0) = 1 , 

(zz) unitary : E*(9) = E*(—9) , 

(in) positive definiteness: Ylkk '=i 'ZkZk'e~ lIm ^ k9k '^ 2 E„(9 k — 9 k ') ^ 0 for any K ^ 1 and 
z k E C (k — 1, 2,..., K) , 

(■ iv ) regularity: the continuity of the map r i —> D(r9) implies that the function r H- h,(r 9) 
is also continuous. 

Then by the Araki-Segal theorem, the properties (i)-(iv) guarantee the existence of the 
unique normal state u> p s over the CCR algebra ^(M’s) such that E*(9) = u> p s(w(9)). 
Taking into account (a) and (j2.16jl we conclude that in contrast to the initial state oj$ the 
limit state oj p s is gauge-invariant. 

(c) The convergence (12.14)1 can be extended by linearity to the algebraic span of the set 

of Weyl operators {iu(a!)} aec . Since it is norm-dense in C*-algebra g/(2rf?s), the weak*- 
convergence of the states to the limit state u> p s follows (see |BRlj . |A,JPlj ). □ 


10 





































Remark 2.4 (a) By Theorem 13.31 (a)-(b), one has pf = pf(r), i.e. the limit state u> p s 
is invariant under the one-step evolution Tf 0 . Comparing (| 2. 7p and \2.16\) one finds that 
pf P*{ v ) f or 0 < v < t. Instead, the evolution for repeated perturbation yields the 
asymptotic periodicity: 




for t = (n — 1 )t + v(t) 


(2.17) 


(b) Let p\ in 12.5\) correspond to the quasi-free gauge-invariant Gibbs state for the inverse 
temperature /3 > 0 and let co Po (-) be any initial normal state of the subsystem S. Since 


u Pl (w(9)) = exp 




holds, we obtain for \ 2 . 14 \) : 


D{9) = exp 


1 |g| 2 

4 1 - \g ry z c \ 2 



(2.18) 


(2.19) 


Put A ct (t) := |g°M; CT | 2 (l — \g a z° r \ 2 ) 1 G [0,1) (Remark II. PI) . Then for the characteristic 
function of the limit state in Theorem \2.tA we get 


^*(^(69) = ex P 



f(l —A"(r))^±^± + A"(r) coth§) 
V a.-cr + 2) 


( 2 . 20 ) 


If w a = 0 (i.e. A <T (r) = 0), the subsystem S seems to interact only with reservoir IZ, 
and it evolves to a steady state with characteristic function 


E* 0 (9) = exp 


|P | 2 <r_ + cr+ 
4 cr_ — <J + 


0 < a + < cr_ , 


( 2 . 21 ) 


which corresponds to the quasi-free Gibbs state for the inverse temperature /3*o : = ln(cr_/cr + ). 
This reflects thermal equilibrium between S and IZ. In this sense, /3*o is the inverse tem¬ 
perature of the external reservoir IZ INVZlj . 

If w a 7 ^ 0, the steady state \2.20 1) of subsystem S has the characteristic function 


E* ( 9 ) = exp 




( 2 . 22 ) 


where the inverse temperature fif (r) is defined by 

coth * ^ - = (1 — A CT (r)) cotli —- + A°Yr) coth — . 
2 v v )) 2 \ ) 2 

Note that filfir) satisfies either fi^ ^ fif (r) ^ fi or /3* 0 ^ fi*( T ) ^ fi¬ 


ll 



































2.2 Correlations: subsystems S + Sn and Sm + Sn 

To study quantum correlations induced by repeated perturbation, we cast the first glance 
on the bipartite subsystems S + S n and S m + S n . We consider the initial density matrix 
( 12 .5p satisfying 


w po(™( 0)) = ex P 


\ 0\ 2 ,, A> 
~ coth Tj 


, uj p Aw(6)) = exp 


l^l 2 P 

~ coth 2 J 


(2.23) 


From (11.201) and (11. 371) . we have: 

Proposition 2.5 For evolved density matrix p(Nr) = Tff TO p the characteristic function 
of the state oj p (nt)(-) is 


^(Nr)(W( 0) = (p | T’ N l„(W(Q)) x = exp 


i(C,X'(Afr)C> 


(2.24) 


where X a (Nr) is the (N + 1) x (N + 1) matrix given by 


X u (Nt) 


ttct * 
v A 


'N 



xUf...Ufj + 


<7- + °+ , 

<7_ — (J _|_ 

<7- + <7± j 
cr_ — cr+ 


1 + e ^ 

1 — e _/? 


/ + 


1 + e - ^ 0 
1 - e~P° 


l + e- p 
1 — e~h 



(2.25) 


Remark 2.6 In the theory of quantum correlation and entanglement for quasi-free states 
the matrix X a (t) is known as the covariant matrix for Gaussian states, see / Ad Ilf . \Ke\J . 
Indeed, differentiating \2.2f | ) with respect to components of ( and ( at ( = 0, one can 
identify the entries of X a {t) with expectations of monomials generated by the creation and 
the annihilation operators involved in A 1.281) . A1.2IA) . 


Subsystem S + S n . For 1 < n ^ N the initial state ^ 5+5 (•) on the Weyl C*-algebra 
srf(J#o) <g) £/of this composed subsystem is given by the partial 
trace 


71 — 1 N 

w 5+sJ^(“o) ® w(«i)) = cj p (w(o 0 ) ® ® 1 ® w(«i) ® (g) 1 ) 

F—i k=n -\-1 


= exp 


N 


fcoth — 
2 


exp 


ai 


k =1 

teeth 2 
2 


(2.26) 

This is the characteristic function of the product state corresponding to two isolated sys¬ 
tems with different temperatures. Put ^ 0 ’ n ) := t (ao, 0,..., 0, «i, 0,..., 0) G C N+1 , where 
ai occupies the (n + l)th position. . Then we get 

"£&.(»(« 0 ) ® = u^ Nr) (W( C ( “’">)) ■ (2.27) 

For the components of the vector Uf ... b r ^C ( ' 0 ’ n \ we get from Remark 11.21 that 

(U? ...Uff C (M )fe = (2.28) 
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'e iNre [(g a z a ) N a 0 + ( g a z c ) n ~ x g a w a «i], (k = 0 ) 

e iNTe [(g a z (T ) N ~ k g a w rT ao + (g a z a ) n - k ~ 1 (g a w C7 ) 2 «i], (1 ^ k < n) 

e tNre [( g a z a ) N ~ n g a w a a 0 + g <T z a (—r) aq], {k = n) 

yiVre (g° z °)N-k g <T w v ao (n<k^N). 

Substitution of these expressions into (I2.24p and (I2.25[) allows to calculate off-diagonal 
entries of the matrix X a (Nr) for £ = which correspond to the cross-terms involving 

«o and op. 

Because of \g <T z a \ < 1 (Remark ll.2j) . these non-zero off-diagonal entries will disappear 
when N —y oo for a fixed n. Hence, in the long-time limit the composed subsystem 5 + iS n 
evolves from the product of two initial equilibrium states (j2.26[) to another product-state. 
On the other hand, the cross-terms will not disappear in the limit N, n —> oo, when N — n 
is fixed m It is interesting that in this case the steady state of the subsystem S keeps 
a correlation with subsystem S n in the long-time limit. 

Subsystem S m + S n . We suppose that 1 ^ m < n ^ N. Then the initial state w(j m+< s n (-) 
on stf <S> — stf (^n) <S> srf (J^ n ) of this composed subsystem is given by the partial 

trace 


771—1 71—1 N 

w 5 m +s„(^( a i) ® w(a 2 )) = 1® w(ai) ® (x) 1 <g> w(a 2 ) <8> (££) 1) 

k =0 k=m +1 k=n +1 

12 Pi r H 2 p 


= exp 


coth 2 
4 2. 


exp 


coth — 


(2.29) 


This is the characteristic function of the product-state corresponding to two isolated 
systems with the same temperature. 

We define the vector *((), 0,..., 0, aq, 0,..., 0, a 2 , 0,..., 0) E C N+1 , where aq 

occupies the (m + l)th position and a 2 occupies the (n + l)th position, then 

C (m -”>)). (2.30) 

Again with help of Remark 11.21 we can calculate the components of f/f... as 


(UZ...UZ C (m ’ n) ) fc = ( 2 . 31 ) 

' e iNre (< g a z a ) m ~ l g a w a [oti + {g a z a ) n ~ m a 2 ] (k = 0 ) 

e iNre (g a z (T ) rn ~ k ~ 1 (g CT w a ) 2 [e*i + (g a z a ) n - m a 2 ] (1 ^ k < m) 

e lNre [g a z cr (—T) aq + ( g a w a ) 2 (g cr z°') n_m_1 a 2 \ (k = m) 

e iNre ^g<r z a^n-k -1 ^ g a w a^2 a ^ (m < k < Tl) 

e iNre g° z cr(_ T ^ q , 2 (/j — n ) 

0 (n < k ^ N) 


The correlation between S m and S n , i.e. the corresponding off-diagonal elements of 
X a (Nr) are non-zero when w ^ 0, and large for small n — m and they decrease to zero 
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as n — m increase. Note that in contrast to the case S + S n (j2.28f) the last components 
n < k ^ N in (j2.31|) as well as the state (I2.30p do not depend on N. This reflects the fact 
that correlation involving S m and S n via subsystem S is switched off after the moment 
t = nr. If w — 0, then (j2.31|) implies that A" CT (Vr) is always diagonal and that dynamics 
(I2.30P keeps S m + S n uncorrelated. 


3 Time Evolution of Subsystems II 


The arguments of Section [2721 indicate that the components in the subsystems S + Sjy + 
... + Sjv-n have large mutual correlations for small n ait — Nr even when N large. And 
those correlation seems asymptotically stable as N —> oo. 

In this section, we consider the correlation among those components simultaneously for 
product initial densities. For this aim, let us divide the total system into two subsystems 
S H) k and C Ut k at the moment t = hr, where 

Sn,k — S + Sk + Sk-l + . . . + <Sfc_ n _|_i, (3-1) 

and 

C n ,k — Sn + ... + <Sfc_(_i + Sk~n + ... + Si. (3.2) 

Here, n € N is supposed to be fixed small and N e N large enough. We may imagine that 
the “cavity” S and “atoms” S\, ... ,Sn are lined as 


5 


N, 


> S : Sk , • • • j Sk—n+li Sk—m • • ■ j S\ 


at this moment. The interaction between S and each of Si,... ,Sk has already ended, and 
they are correlated. While Sk+i, ■ ■ ■ ,Sjy have not interacted with S, yet. Let us regard 
that S nt k is the “state” at t = kr of the time developing single object S^ n . That is, S^ n 
has S,Sk, ■ ■ ■, Sk—n as its components at the time t = kr. And it develops changing its 
components as well as the correlation among them. As the time pass from t — {k — l)r 
to kr , the “atom” Sk enters into S^ n and the “atom” Sk~ n leaves from S^ n . It is also 
possible to regard S^ n as the view from the window which is made to look the “cavity” 
and the n “atoms” just have interacted with the “cavity”. 

We are interested in S^ n , since it might be interpreted as a simplified mathematical 
model of physical objects in equilibrium with the reservoir or of metabolizing life forms 
which maintain their life by interacting with the environment, i.e., the macroscopic many 
body systems which are macroscopically stable but exchange their constituent particles 
as well as energy with the reservoir microscopically. 

Below we consider the large-time asymptotic behavior of state for S^ n , i.e., for the 
subsystem S ri) k with fixed n and large and variable k for the initial state (12. 5 p with general 
density matrices po,Pi £ <£ 1 ( 2 ^). 

To express the state of S^ n at t = kr, we decompose the Hilbert space Jif into a 
tensor product of two Hilbert spaces 

^ = ^s n , k <S)^c n , k . 
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Here Ms n k is the Hilbert space for the subsystem (13.11) and d%?c nk f° r (13.21) : 


k 



j=k-n -\-1 



(3.3) 


If p G is the initial density matrix of the total system S n ^ + C n ^, the reduced 

density matrix ps^ n {kr) of at t = kr is given by the partial trace 


Ps~ n (kr) = Tijp Cnk {T^ Q p) = Tr^ ci (Tr^(T^p)) , 
for in (12.2)1 . where we decompose Mc n k as 


k—n N 

^ (g) ^c 2 , ^ = (g) ^ , ^c 2 = ® 

J = 1 i=fc+l 


(3.4) 


3.1 Preliminaries 

Here we introduce notations and definitions to study evolution of subsystems in somewhat 
more general setting than in the previous sections. 

In order to avoid the confusion caused by the fact that every M’j coincides with & 
in our case, we treat the Weyl algebra on the subsystem and the corresponding reduced 
density matrix of p G in the following way. On the Fock space J r ®( m + 1 ) for 

m = 0, 1 ,..., N, we define the Weyl operators 

W m ( C) := exp (i (C,b)m + i-HbX)m+i ^ ( 3 . 5) 

where ( G C m+1 , bo,... , b m and Sq, ..., b* m are the annihilation and the creation operators 
in J r ®( m + 1 ) ) which are constructed as in (11.31) satisfying the corresponding OCR and 


(C j b) m+l ^ ^ C jbj i 

3=0 


(b,O m+ i = J2&'r 

3=0 


By we denote the C*-algebra generated by the Weyl operators (13.51) . 

To discuss the dynamics of our open system, it is convenient to introduce the modified 
Weyl operators (cf. Proposition 11 .11) 


W"(C) := exp 


' cr_ + G + 

-4(u_ - o+) 


(C) C)m+1 


W m (o , 


(3.6) 


for m — 0,1, 2, ..., where ( G c m+1 and (•, • ) n+ i denotes the inner product on C m+1 . We 
also use the notation 


w a (6) := W°(6) for 6 G C. 


(3.7) 
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Below, we adopt the abreviations: 


<c/ (m) = and %? (m) = £!(i^ (m+1) ) (3.8) 

for the Weyl C* algebra on t ^'®( m + 1 ) and the algebra of all trace class operators on &®( m+1 ) 
for m — 0,1, 2, ..respectively. Note that the bilinear form 

( • | • ) m ■ tf (m) X ^ (m) 3 (p, A) hg Tr [pA] G C (3.9) 

yields the dual pair Indeed, the following properties hold: 

(i) (p | A) m = 0 for every A G implies p = 0; 

(ii) ( p | A) m = 0 for every p G ^ m ) implies A = 0; 

(iii) |(p|^)m| < IMkll^lU- 

These properties are a direct consequence of the fact that is weakly dense in 

£(j^®( m +!)) the dual space of Below we shall use the topology in¬ 
duced by the dual pair (^ m ), on ^ m ). We refer to it as the weak *-g/ ( ' m '> topology, 

see e.g. |Ro] . pBRlj . 

Note that for the initial normal product state (12.51) the calculation of the partial trace 
over in (13.41) is straightforward: 

N k 

Tr* 2 (Tpy ® ft) = (g) ft . (3.10) 

j =o j =o 

Here stands for the evolution map (11.181) on for k ^ m ^ N. 

To check (13.101) . it is enough to show 

N k 

® ftimK)® i>« = (itf ®Pj|w*(o>* (3.ii) 

3=0 j =0 

for any ( G C fc+1 , where 1 is the unit in algebra g/( N ~ k ~A. Let ( G C JV+1 be defined by 
Cj = Cj f° r 0 < j < k, Q — 0 for k < j ^ N. Then W k (() ® 1 = Wjv(C) holds. Remark 
11.21 readily yields 

ui [N) ... u a k {N) c = (ui {k) ... ui {k) o~ . 

Together with (II,37jl . it follows that 

I’m *‘( ^A-(C)) = (c!? - into) ® i. 

which implies 


N 


((g)/b 

J=0 


rpCJ ( N )* 
/CT ,0 


(W*(C) »!)}« = 


(® ft I i£S )- w*(c»*. 


J=0 


(3.12) 
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This proves (13. lip and thereby the assertion (13.10j> . 

Here we have used the notation U a ( } k) for the (fc + 1) x (k + 1) matrix whose components 
are given by 


e %Tt g a (t){ 5joz a (t) + 5 je w a (r)) (■ i = 0) 

( u e (k) )ij = { e ire g a (T){S j0 w a (T) + 6 je z a (-T)) (i = £) 


(3.13) 


e iT % 


(otherwise) 


for £ = 1,2,... ,k (c.f. Remark fl.2l ). Then the one step evolution T[ (k> on is given 
by 

(T/ (fc) p | W k (0)k = (P I T^ k) *W k {Q)k 

where 


T^ k) *W k ( C) = exp 


er_ + cr + 


■«C,C)fc+i - (U? {k) (,U^ m \) k+ 1 )]w k (U^ k \) , (3.14) 


4(a_ - a+) 

p G and ( G C fc+1 (see Proposition II.ip . 

To calculate the partial trace (13.4p with respect to J^ C1 , we introduce the imbedding: 

/Co \ 


r m+l,m ■ C m+1 9 C = 


/Co\ 

Cl 

C 2 


\C m) 


0 

Cl 

C 2 


\C m) 


— l"m+l,mC ^ C 


m+2 


(3.15) 


for m — 0,1, 2 ,..., N and the partial trace over the second component R rn .m.+i '■ ->■ 

^ m ) characterised by 


(iC,m+lP|w(Co)®w(Cl)<S>-- ■®w{Cm))m = (p|t/>(Co) ® 1 ® ttf(Cl) ® ® w(Cm))m+l (3.16) 

for p G ^ m+1) . Therefore, its dual operator R* n m+1 has the expression: 

^+1^(0 = ^m+i(r m+ i, m C) for C e C m+1 . (3.17) 


Lemma 3.1 For m G N and £ = 1,2,... ,m, 

U°i7 +1) r m+ i, m = , (3.18) 

holds. 
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Proof : In fact, for the vector ( = *(Co, Ci,•'' , Cm) € C m+1 , one obtains 

(U!iT'’rm+l, m Qi = (r„, +h „,UZ (m) Oj 


e-V(r)(^(r)Co + ^(r)0) 

U = o) 

0 

(J = 1) 

e‘%-1 

(2 ^ J < £) 

e iTe g a {T) (w a {r) Co + z c {-t)Q) 

U = * + 1) 

e^0_i 

(£ + 2 ^ j ^ m + 1) 


by explicit calculations. This proves the claim (13.1 8[) . □ 

For k E N and m = 0,1, 2,..., k — 1, let the maps rk, m ■ C m+1 —> C fc+1 and R m ,k '■ 

<g{k) 

—> be dehned by composition of the one-step maps (13.1511 . (13.1611 : 


and 


^k,m ^k,k—1 ° ^k—l,k—2 ° • • • ° ? 


Rm,k Rm,m +1 ° Rm+l,m+2 ° • • • ° Rk—l,k ? 
respectively. This definition together with ( 13.161 ) and (13.171 ) imply that R^ k 
and 


g/W 


R m,k £(Co) ® ty(Cl) ® ® ty(Cm) = w(Co) ® 1 ® . . . <8> 1 <8> w(Cl) ® ® w(Cm) • (3.19) 

Hence, by ()3. 16(1 the map R m ,k, which is predual to (13.1911 . acts as the partial trace over 
the components with indices j = 1, 2,..., k — m of the tensor product <S> k =0 Pj E 
Therefore, the map R, h k coincides with the partial trace Tr Cl in (I3.4jl . Then R 7l) k combined 
with (13.1 Op gives the expression 

k 

Ps~„(kr) = Rn,kT^((^) pj ) for k ^ n + 1. (3.20) 

3=0 

We summarise the action of the above maps (13.1511 . (13. 1711 - (13. 1411 on the modified Weyl 
operators (13.611 by 

Lemma 3.2 Let k E N. Then, 


(i) 

R m,m+1 (y^m (C) ) ~ ^m+l ( r m+l,mC) j 

(3.21) 

(ii) 

-^m,m+fc(^m(C) ) = m+fc( r m+fc,mC) 

(3.22) 

holds for m 

= 0,1,2,..., C e c m+1 ; and 


(hi) 

r; (m) *(w , "(C)) = , 

(3.23) 

(iv) 

m). 0) = ... c/“i7 ) c/“ (m) o, 

(3.24) 

(v) 

Tjcr(m+k) _ jj<y{rn) 

^l-L-k 1 m+k,m ' m-\-k,m {J £ 

(3.25) 


holds for m E N , C £ C m+1 and £ = 1, 2, ..., m. 
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Note that the claim (v) in the above is an obvious extension of Lemma 13.11 This lemma 
yields the following statement. 

Lemma 3.3 For m,k E N ? l = 1,2,... ,m, 




7 ~>* /yiCr(?77.-|-fc)* D* 

^m,m+k-*- Z ^ £+k ^m,m-\-k 

(3.26) 

holds 

on s/t™) _ Therefore 

t~> ry-icr(m+fc) _ rpCr(m) p 

£-|-/c ^ ff ^m,m-\-k 

(3.27) 

and 


p rp(r(m+k)rpcr(m+k) rpcr(m+k) 

-tlrriim+k- 1 - k 1 k-1 •••■*■ l 



_ ( 73 rpcr(m+l) \ 

— (rXm,m+ 1^1 ) 

/ T-) r 7 -icr(ra+/c — 1)\ / 7~> rr-,cr(m-|-/c) \ 

• • • v-^m+fc— 2,m-\-k—1 -*-1 / l,m+fc-* 1 / 

(3.28) 


hold on ^h m+fc ). 


Proof : The identity (13.261) follows from Lemma 13.21 by considering the action on the 
modified Weyl operators. By taking its adjoint, (I3.27|) follows. A simple application of 
induction over k yields the last identity. □ 

Let us concentrate on the evolution of the subsystem S, first. To this aim, we introduce 
the map T[ ■ \ ■} : x —>■ to express the one-step evolution 

T[p 0 \p 1 ]=RonT^ 1 \p 0 ®p 1 ) for p 0 , Pl e^ (0) , (3.29) 

of the density matrix p 0 under the influence of p 1 , see (13. 201) . We also denote by 

T\p\ ■= e -^a*a pe iera*a for p £ tf(0) , ( 3 . 30 ) 


the “/ree” one-step evolution of density matrix corresponding to any of subsystems 5*,, 
c.f. (jl.8[) . Then one obtains the following assertion. 

Lemma 3.4 For any l, m E N fulfilling £ ^ m, £ E C m , 9 E C and p 0 , pi,..., p e E 
the following properties hold: 

r*®m K _ i(c) ] = w^fie^C), (r- 1 )*® m [W^_ 1 (C)] = W^fie-^C) ; 

TRo,i = AoyT® 2 ; 

(T[po\pi] I w a {9 ))o = (p 0 | w C7 (e ieT g a (T)z a (T)9)) 0 (p 1 | w a (e ieT g a (T)w a (T)9 )) 0 ; 

r(r[p„M) = t{t Po \t P i ], r-yTpoipj) = r[r-Voir-Vi]; 

R,- U T;‘"[p a ® pi ® ... ® /)f] = T[p 0 |pi] ® Tpj] ® ■ ■ ■ ® T[p.j. 

//ere (.)®( m+1 ) denotes the (m + 1 )-fold tensor product of the corresponding operator (•). 


(in) 

(iv) 

(v) 

(vi) 
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Proof: (i) Since 


<p|r*K( 0 )]}o = {Tp|K( 0 )]}o 

= Tr(pe“™‘“S' , (»)e-" T “' a ) = H:(pffi"(e“ T 0)) = (p | {S"(e‘ cr g)) 0 

holds for p G we obtain the desired equality for m = 1. The equalities for m > 1 
follow from (11.271) and the definition of tensor product of 

(ii) Taking into account (i) of Lemma [3.21 and the above (i), we get 

( TR 0 , 1 )*w (7 {9) = R* 01 T*w a (9) = R* 01 w"(e ieT 9 ) = W{ (e ieT r lfi 9) 

= T m *W?(r lfi 6) = T m *Rl 1 w C7 (9) = (R 0jl r^ 2 )*w (T (9). 

(iii) By virtue of (13.151) and (j3.13[) one obtains 


Ul {1) r lfi 9 


/ e ire g a z a^ e ire g a w a \ / q\ 

\e ir( g rT w rJ , e iTe g a z a (-r) J 


/ e ire g a z aQ\ 

~ y e ire g a w aQJ , 

which impllies 

(T[p 0 |pi] | w"(9 ))o = (i? 0 ,iTf (1) [p 0 ® Pl ] | w*(0)) o 
= (Po ® Pi I 7T (1) *^, 1 ^(^)) 1 = <Po ® Pi I bhf (7/f (1) r li0 P))i 

= (p 0 | w (7 (e ieT g a z <T 9))o(pi \ w"(e^g"w°9)) 0 . (3.31) 

(iv) By applying the adjoint operators 


to the modified Weyl operators, we get 

|''y _ ±l^(g>(r7i+l)*2~ ,<T ( m )* _ j-icr(m)* ^<y-±l^(g)(?ra+l)* 

from (iii) of Lemma [3.21 and the above (i). Then, duality derives the assertion. 

(v) These identities follow from the above (ii), (iv) and the definition (I3.29p . 

(vi) Let ( G C e . By (13.15jl and (13.131) . we obtain 

tff ( %-iC = e iT£ \g° f Co, 9°w a Co, Ci, ■ • • , Ct- 1), 
where *(•••) is the vector transposition. Then we get 

{Re-i,eT^ e \p 0 <S> pi ® ... ® p^) | kb/L^O^-i 

= (po ® pi ® ■ ■ ■ ® pt I rf w *i?|_ v w7_i(c))*! = (p 0 ® pi ® ■ ■ ■ ® p< I w7(c/f w r^_iC))^ 

= (Po | {e leT g a {r)z a { t)Co))o{pi \ w a {e leT p <7 (r)w ,7 (r)Co))o(p 2 | w CT (e* eT Ci))o • ■ ■ (fit I w ff (e* er O-i))o 
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= (7>o|p.l I ffi"(Co))o{p 2 1 T‘ K(Ci)]> 0 ... (ft I r K(Cfti)]>. 

= (TVqIpi] ® T[p 2 ] ® ■ ■ ■ ® T[ft] | W£. i(C))e-i, 

where we used (iii) and (i) at the fourth equality. These finish the proof of the lemma. □ 
Note that (13.311) coincides with (12.71) . 

Next, we consider the multi-step evolution for the subsystem S. To this aim, we define 

7"(*0 ; <^(0)(fc+l) <^(0) jgy 

T (fc) [Po|pi, ...,Pk] = Ro,kT^ k) T^ ... Tf {k \p 0 ® Pl <g> ... <g> Pfc ), (3.32) 

for k G N and p 0 , p 1 ,..., p k G < ^’^, c.f. (|3.20j) . (13.291) . 

The following lemma holds. 

Lemma 3.5 For any 6 G C, fc G N and p 0) Pu ■ ■ ■ •> Pki ■ ■ ■ G < ^°\ £he following properties 
hold: 

(i) T<‘>[ft|ft] = r(po|ft] ; 

(ii) r< i+1 )[p„|p 1 .....p t+ i] =r[r(‘)[p„|ft,...,p*]|rWi] 

= r w [r(ft,|pi]|rp2,....rp t+ i] ; 

(iii) r(r«(po|Pi,....Pd) =TM[Tp„|Tp,,...,Tft] ; 
r- 1 (r<‘ ) [p„|ft....,p l ]) =n‘>[r-y„|r-y 1 ....,r-y i ] ; 

(iv) r t ’ ( ‘ +m) ... rr l ' I+ ’” ) |po ® pi ® • • • ® p t+m ] 

= T (t| [pci|pi, ■ ■ ■ ,pt] ® T‘[p t+ i] ® ... ® 7*[p*+m] for m = 0,1,2,...; 

(v) (T w [po|pi I .... pj | fiT(0))o = (Po I ff’(e i “ T ( 9 V-")‘0))„ 

k 

x n<ft I «i’(.e ,k ~( g ’zn k -’g a W' 0 ))„. 

3 =1 


Proof: (i) This is obvious by definition, 

(ii) By Lemma [3731 we get 


p rriCr(k^ rTl&(k') rr\&(k') 


l T3 r p <J P)\f TD r n<jyz.)\ /p rpUyKj\ 

f-Ko.lTl )\^ 1,2^1 J • • • 1 J 


-,cr(2) 


°-( fc b 


Then, definition (13.32ft and Lemma [T4l' vi) yield 


r ( ‘>[polPi, ■ ■ ■ ,p»] = (flo,irr (1) )(fli,2rr (2) ) • • • (Rk-i,kTf {t) )(p 0 ® Pi ®... ® p t ) 

= Tin •• 7 Tr[p„|p 1 ]|rp 2 ]... |T‘- 2 p l _ 1 ]|T l - 1 p t ], 


which iplies the claim. 

(iii) This can be derived by induction using above (ii) and Lemma 13. 4f iv). 

(iv) Due to Lemma [3731 we have 


R 




-icr(/c+ra) rj-icr(h-\-m) 


k—1 


... t; 


i cr(fc+m) 


(Po ® Pi 


Pk+m ) 
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_ ( TD rricr(m+ 1)\ / t-> rpcr{rn+ 2)\ 

*-771,771+1 -L 1 J ^-^771+1,771+21 / ‘ ‘ ‘ 

... (i? m+ j fc _i )TO+ j t T 1 <T(m+fc) ) (p 0 ® Pi ® • • • < 8 > Pfc +m ). 


Using successively Lemma l3.4T vi) and the result (ii) of the present lemma, one obtains 
the assertion. 

(v) By virtue of Lemma I3.4lf ii) and of the result (i) above, one can prove the case k — 1. 
Let us assume the validity for k ^ 1. Then the validity of the case k + 1 follows from the 
(ii) above and the formula 


<T (fc+1) [p 0 |pi, • • • ,p*+i] | w a (e))o = (T[r (fc) [p 0 |pi,.. .,p k }\r k p k+1 ] I n+0))o 


= <r (fc) [p 0 |pi, ...,Pk}\ w <7 (e ieT g a z a 9))o ( T k p k+1 \ w"(e^ g " w "9)) 0 
= (po\w°(e ikeT (g' T z") k e ieT g"z°0)) o 


x Y[(Pj 

3 = 1 


-^'(7 

w 


(e ikeT (g a z a ) k ~^ g a w r7 e ieT g a z c 


0))o {p. 


'k +1 


I ^(7 

w 


y(" + 1 )eT 9^^)) 0 , 


which proves the assertion (v) by induction. □ 

Here, we comment that Lemma 13.51 (v) is a revisit to the evolution of the subsystem 
S in Lemma [2.21 


3.2 Reduced density matrices of finite subsystems 

In this section, we consider evolution of subsystems S n ^ (13.11) and S~ n . Our aim is to 
study the large-time asymptotic behaviour of their states, when initial density matrix is 
given by (12.51) . 

For the density matrix pi in (12.51) . we assume the condition: 

OO 

[H] D{6) = J^[(pi | w((g a z a ) l 9)) 0 converge for any 9 G C 

1=0 

and the map R 3 t K > D(t9 ) G C is continuous. 

Here, we do not assume gauge invariance of p\. (c.f. Theorem 12.31) 

Under the condition [H], one obtains the following theorem: 


Theorem 3.6 There exists a unique density matrix p* on TF such that T[p* | pi] = Tp* 
holds. And p* also satisfies 

\9\ 2 a_ + o + , |+(rK(r )| 2 


(1) u p fiw(9)) = exp 


1 - 


1 - \g a (r)z a (r)\ 2 


4 o\_ — cr + 

(2) T {k) [p*\p 1 ,...,p 1 ] = T k p*fork>l; 

(3) For any density matrix p 0 in \2.5\) . the convergence lim^oo 
p* holds in the weak*topology on 


D(g*(T)w°(T)9); 

r~ k 


0 | Pi, ■ ■ ■ 5 Pi 
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Remark 3.7 (a) The weak*-srf^ topology on induced by the pair A3.91) 

is coarser than the weak*-C(3P) topology, which coincides with the weak and the norm 
topologies on the set of normal states [Roi \BR1V . 

(b) When p 1 is gauge-invariant, the characteristic function in (1) coincides with A2.1&) 
and the present theorem reduces to Theorem \2.3[ Especially, the free evolution T[p * | pi] = 
T[p*] reduces to the invariance T[p* | pf\ = p*. 


Proof: First, we note that hm^ 00 (po|tc CT ((g cr A T ) fc 6 l ))o = 1 because of \g a z a \ < 1 and of 
the weak continuity of the state u> po = (p 0 | • )o- Then by Lemma [3751 (iii),(v) and Lemma 
I3.4f i). we get 

lim {T~ k (r (t) [p 0 1 Pi,.... Pi]) K(»))o 

= lim (T~ k [po]\w a (e iker (g CT z C7 ) k 0)) o 

/c—>-oo 
k 

X (3.33) 

j=0 


=n^i^(^)y^)>0 

1=0 


= exp 


-\() 2 a + rr. \g a w a \ 2 
.4 cr_ — cr + 1 — |g 0 A CT | 2 


D(g a w a 9) 


5 


which means that lim (T k ( T (k> [po \ p\,..., pi]) | , u; cr (6 l ))o is equal to the right-hand side 

of (1) in the theorem. (Recall (13.6j) and (13.71) .4 The right-hand side of (1) satisfies: (i) 
normalization, (ii) unitarity and (iii) positivity, and (vi) regularity, since it is a limit of 
characteristic functions, under condition [H]. Hence from the Araki-Segal theorem as in 
Section [All there exists a state u;* on the CCR-algebra such that its characteristic 

function is given by the right-hand side of (1). Moreover, the continuity assumption about 
the function D yields that the state cu* is normal by the Stone-von Neumann uniqueness 
theorem |BR2j . Hence, there exists a density matrix p* such that u;* = a; Pt , which conclude 
(1). Now, (3) is obvious. 

Free evolution T[p* | pi] = Tp* can be derived from (1) by the use of Lemma 13.41 
(iii),(i) and (13. 6 1) . (13.7|) . Indeed, one has 


{T[p*\pi}\w a (9)) 0 = (p*\w a (e ieT g a z a 9)) 0 ( P i\w a (e i€T g a w a 9)) 0 
O- + er+ 


= exp 


(\g°z"9\ 2 + \g°w°9\ 2 ) 


= exp 


4(a_ - o + ) 
x (p*\w(e ieT g a z a 9)) 0 (pi\ui(e i<:T g a w a 9)) 0 
cr- + a + , IP^I 2 , _ ct „ |2 


(3.34) 


g a w a 9\‘ 


4(cr_ — CT + ) V 1 — \g a z a \ 2 

xD(g a w a e ieT g CT z a 9)(p 1 \w(e ieT g a w a 9)} 0 


\g a z a 9[ 
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= exp 


er_ 


= exp 
-0 + 


cr_ 




<j n ..cr 12 




-D(e ieT g a w a 9) 


.4(cr_ — cr_|_) 1 — |g 0 'z°'| 2 
|0| 2 1 <p*|£(e^)> 0 = <T[p*]K(h)> 0 = uW^(*)) , 


L4(a_ - o+) 

where we used the equality D(g a z a 9)(pi\w(9)) 0 = D{9). 

Now the assertion (2) follows directly from T[p* | pi] = Tp*, by using Lemma l331f i0 (iii). 
To prove the uniqueness of p*, let p* be another density matrix satisfying T[p* | pi] = 
Tpt- Then, p* satishes the property (2) and 

p* = lim T~ k [T [k) [p 0 | Pi, - - -, Pi]] , 

k—> oo 


which coincides with p* by (3). Hence, one gets p* = p*. □ 

Now we consider the large-time behaviour of the states (13.41) of subsystems S~ n . Let 
pi be a density matrix on & satisfying the condition [H], Then we have the following 
theorem. 


Theorem 3.8 For any density matrix po on FF and n,m E N, m ^ n, the limit: 
(r -* ) ® (m+ i) flmm+ ^«) ) ( po 0 pfl"**)) T$$(p, ® pf”) as oo , 

holds in the weaktopology on < &( m ). Here p* is the density matrix on FF given in 
Theorem 13. hi 

Proof: By Lemma [3.31 Lemma I3.5f iv) and Lemma I3.4f ivl. we obtain 

('T-k\®{m+l) jd rpcr(m+k) / ®(m+fc)\ 

V' J *Fm,m+k-L (n+k)T,0 VP° ® Pi J 

= (T- fe )® (m+ i)T^ m )... Tf (m) i? mim+fe T* (m+fe) ... T? (m+k) (p 0 ® pf (m+fc) ) 

= c T~ k )® {m+1) T ... T^ m) (r (fc) [p 0 |pi,..., pi] ® (r fc [pi]) 0m ) 

= TZ {m) ... T? m \T- k )« m+1 > (T (fc) [p 0 |pi,...,pi] < 8 > (T fc [pi])® m ) 

= (T~ k [T ik) [p 0 \p u ... ,pi]] ® p? m ). 

Since one has 

hmr-‘(r<‘>[p„|p 1 ,...,p 1 ])=p. 

fc—»0O \ / 

in the weak*-.^ 0 ) topology, we obtain also the weakconvergence 

(r- k [T^[p 0 | PI, • • •, Pi]]) ( 8 ) pf m —► p* ( 8 ) pf m as fc oo . 

By the duality (13.91) . one also gets the continuity of T^'g* and hence, the weak*-^/^ m ) 
convergence 

2^?(7^‘[T ( ‘>[po|pi,....Pi]]®pf”) ^ Tffl (p, ® p?™) ask^oo, 

claimed in the theorem. □ 

Let us put m = n in the theorem. Then by (13.201) . we obtain the limit of the reduced 
density matrix ps~ n {-) for the subsystem S~ n : 
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Corollary 3.9 The convergence 

lim (r- t )® ( “ + 1 V 5 ™((« + k)r) = T^(p, ® pf) (3.35) 

rC—^OO 

holds in the weaktopology on 

Since T is the free evolution (13.301) , the limit (13. 3 Of) means that dynamics of subsystem 
S^ n is the asymptotically-free evolution of the state, which is given by the n-step evolution 
of the initial density matrix p* <g) pf n of the system S + C n . 

From the continuous time point of view, the subsystem S^ n shows the asymptotic 
behaviour, which is a combination of the free evolution and the periodic evolution, c.f. 
Remark 12.4( a) . 
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